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It is demonstrated that strong laser pulses can introduce torsional motion in the axially chi- 
ral molecule 3,5-dinouro-3',5'-dibromo-biphenyl (DFDBrBPh). A nanosecond laser pulse spatially 
aligns the stereogenic carbon-carbon (C-C) bond axis allowing a perpendicularly polarized, intense 
femtosecond pulse to initiate torsional motion accompanied by a rotation about the fixed axis. We 
monitor the induced motion by femtosecond time-resolved Coulomb explosion imaging. Our theo- 
retical analysis corroborates the experimental findings and on the basis of these results we discuss 
future applications of laser induced torsion, viz., time-resolved studies of de-racemization and laser 
controlled molecular junctions based on molecules with torsion. 



I. INTRODUCTION 

Gaining control over the external and the internal de- 
grees of freedom of molecules is of great interest for a 
range of areas in molecular science. One approach, dat- 
ing as far back as the 1920s [TJ [2], relies on the applica- 
tion of inhomogenous (static or low frequency) magnetic 
and electric fields. This has proven most useful for con- 
trolling the full three-dimensional motion of molecules, 
including control of the velocity distribution, bringing 
molecules to a stand still, storing them for extended pe- 
riods of times and trapping them in a confined volume 
(for a recent review see [3]). An alternative and comple- 
mentary approach relies on the use of non-resonant non- 
ionizing laser fields typically supplied by pulsed lasers. 
Many studies during the last decade have shown that 
strong non-resonant laser fields can effectively manipu- 
late the external degrees of freedom of isolated gas phase 
molecules. The manipulation results from laser- induced 
forces and torques due to the interaction between the in- 
duced dipole moment and the laser field itself. Examples 
of manipulation include deflection [4], focusing [5 and 
slowing [6] of molecules through the dependence of the 
non-resonant polarizability interaction on the intensity 
distribution in a laser focus. Likewise, the dependence of 
the induced dipole interaction on molecular orientation 
has proven highly useful for controlling the rotation of a 
variety of molecules [7]. In particular, the spatial orien- 
tation of molecules can be sharply confined with respect 
to axes that are fixed in the laboratory. 

Molecular manipulation by induced dipole forces ex- 
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tends beyond the external degrees of freedom and has 
also been demonstrated for the internal degrees of free- 
dom such as vibrational motion [8] in molecular hydro- 
gen. Notably, the electrical field from laser pulses can 
modify energy potential barriers such that photoinduced 
bond breakage of a small linear molecule is guided to 
yield a desired final product [9j [10] . 



In the case of larger molecules control of the lowest fre- 
quency vibrational modes attracts special interest since 
some of these modes correspond to motion along well- 
defined reaction coordinates separating two conforma- 
tional minima (conformers). Although many molecules 
contain large number of conformers it is often just two 
conformers that dominate important chemical properties, 
for instance chirality. A particularly important exam- 
ple is found in axially chiral molecules [TTJ [12] such as 
biaryl systems. In these molecules rotation about a sin- 
gle stereogenic carbon-carbon (C-C) bond axis changes 
the molecule from one enantiomer into the opposite enan- 
tiomer (mirror image). Recently, we demonstrated that 
the laser-induced nonresonant polarizability interaction 
can also be used to influence the internal rotation of an 
axially chiral molecule around the stereogenic C-C bond 
axis [13]. In particular, we showed that by fixing the C-C 
bond axis of a substituted biphenyl molecule in space, us- 
ing laser induced alignment by a long laser pulse, it was 
possible to initiate torsional motion of the two phenyl 
rings by a short laser pulse polarized perpendicular to 
the fixed axis. The purpose of the present work is to ex- 
tend our recent paper [13]. In particular, we discuss our 
theoretical modelling in detail. 



2 



II. OUTLINE OF THE STRATEGY FOR LASER 
CONTROLLED TORSION 

A ns laser pulse aligns [7 the stereogenic axis of 3,5- 
diflouro-3\5'-dibromo-biphenyl (DFDBrBPh) along its 
polarization. This enables a much shorter fs pulse, which 
we will refer to as the kick pulse, to both initiate tor- 
sional motion and set the molecule into rotation around 
the fixed axis. The internal as well as external rotational 
motion is monitored by fs time-resolved Coulomb explo- 
sion imaging. 

The inset of Fig. [I] shows a model of the DFDBrBPh 
molecule with the stereogenic axis marked by red. Ac- 
cording to our quantum chemical calculations, the laser- 
free torsional potential, illustrated by the red dotted 
curve in Fig. [I] has minima at dihedral angles of {<pd) = 
±39° corresponding to the R a and S a enantiomeric forms, 
and where the dihedral angle is the angle between the two 
phenyl rings. The twisted equilibrium shape is charac- 
teristic of biphenyl compounds. The traditional view is 
that this non-planarity results from competition between 
stabilization via conjugation of the 7r-orbital systems and 
steric repulsion between the ortho-positioned atoms [14] . 
We might mention that alternative explanations have re- 
cently been discussed [TSJ [I6j [17] . Our strategy for con- 
trolling (j)^ relies on a transient modification of the field- 
free potential curve by a fs kick pulse (Fig. [T]). The mod- 
ification, caused by the nonresonant polarizability inter- 
action [7J [9] [18] converts the initial stationary quantum 
states, localized near the minima of the torsional poten- 
tial, into vibrational wave packets, i.e., coherent super- 
positions of several quantum states. The temporal evo- 
lution of these wave packets gives rise to time-dependent 
torsional motion, and the grey and the black curves in 
Fig. [I] show the calculated expectation values of (j)d for a 
particular set of laser parameters. The stereogenic C-C 
bond axis is fixed in the laboratory frame by adiabatic 
alignment utilizing a moderately intense, linearly polar- 
ized ns laser pulse [7] [19] . The ns pulse is intense enough 
to keep the C-C axis tightly confined, yet weak enough to 
modify the torsional potential only slightly. The fs kick 
pulse is applied with its polarization perpendicular to the 
aligned C-C bond axis to ensure primarily influence on 
torsion while avoiding other possible vibrational motion. 

The DFDBrBPh compound is prepared as a racemate 
and the experiment is carried out on isolated molecules 
at rotational temperatures of a few Kelvin. Under these 
conditions we have equal numbers of molecules initially 
localized in the —39° or 39° conformation (see Fig. [TJ 
and no thermally induced transitions between the two 
conformations occur. The intense probe pulse, sent at 
time t p with respect to the kick pulse, removes several 
electrons from the molecules, thereby triggering Coulomb 
explosion into ionic fragments. In particular, the Br + 
and F + fragment ions recoil in the planes defined by the 
Br- and F-phenyl rings. By recording the velocities of 
both ion species with two-dimensional ion imaging [19] , 
we determine the instantaneous orientation of each of the 
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FIG. 1: (Color online) To illustrate the principle of laser- 
induced torsion and time-resolved de-racemization we show 
the calculated kick pulse induced time-dependent torsional 
potential as a function of time measured with respect to the 
center of the pulse and dihedral angle 4>d between the Br- 
and F-phenyl planes. The asymmetry in the potential is ob- 
tained by orienting the molecules (here with the Br-phenyl 
plane out of the paper), 3D aligning them, and by polarizing 
the kick pulse at an angle of 13° with respect to the second 
most polarizable axis (SMPA) (see Fig. [2}. The red dotted 
curve at large times illustrates the laser-free time-independent 
torsional potential. For the S a enantiomer, starting out with 
{(/>d)i = —39°, the time varying potential induces an oscilla- 
tory motion (grey curve) corresponding to torsion confined 
within the initial well. By contrast, due to the induced asym- 
metry between the two wells, the initial R a enantiomer is 
traversing the central torsional barrier, and ends up as an S a 
enantiomer undergoing internal rotation (black curve). The 
kick pulse has an intensity of 1.2 x 10 13 W/cm 2 and a duration 
of 1.0 ps (FWHM). The torsional motion may be monitored 
by fs time-resolved Coulomb explosion imaging. The inset 
shows a model of the DFDBrBPh molecule with the stero- 
genic axis marked by red (grey). 

two phenyl planes at the time of the probe pulse. 

III. THEORETICAL METHODS 

A. Modelling the dynamics of the molecule 

To account for the motion of the nuclei of the DFD- 
BrBPh molecule we make use of two sets of coordinate 
systems (see Fig. A molecular fixed (MF) frame at- 
tached to the molecule and a laboratory fixed (LF) frame 
specified by the lasers. The MF coordinates are chosen 
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FIG. 2: Model of DFDBrBPh along with the molecular fixed 
(MF) xyz and the laboratory fixed XYZ coordinate axes. 
The dotted line indicates the second most polarizable axis 
(SMPA), which is located 11° from the Br-phenyl ring and 
28° from the F-phenyl ring. 



with the z axis along the stereogenic axis pointing from 
the phenyl ring with the bromines towards the phenyl 
ring with the flourines, and the x axis is chosen along the 
phenyl ring with the bromines. The LF coordinates are 
chosen with the Z axis along the polarization direction 
of the ns pulse and the X axis along the kick pulse po- 
larization direction. In agreement with the experimental 
observations (cf. Sec.[V|), we assume that the stereogenic 
axis of the DFDBrBPh molecule is aligned along the Z 
axis. Also, we neglect all normal modes, except the lo w- 
est one, which corresponds to torsion (see Sec. IIIB for 
details on quantum chemical calculations, justifying, e.g., 
this assumption). In this simplified situation and in the 
absence of the kick pulse the task is reduced to describ- 
ing the coupled rotations of the two (rigid) phenyl rings 
of the molecule as given in the LF frame by the Hamil- 
tonian [Atomic units, m e = e = ciq = /i = 1, are used 
throughout unless indicated otherwise] 
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Here I = Ibt + is the total moment of inertia of the 
molecule for rotation around the stereogenic axis and 
I re \ = IbtIf/I is a relative moment of inertia for the 
two phenyl rings around the axis. A full rotation of ei- 
ther phenyl ring leaves us with the same molecule im- 
plying 27r-periodic boundary conditions of the eigenfunc- 
tions of H mo \ from Eq. 0, i.e., ^(0Br + 27rm, 0p + 27rn) = 
ip(4>Br, 0f)> with m and n integers. We shall assume that 
this property also applies to and </>d, so that we simply 
need to consider eigenfunctions <fid) — £($)x(0d) of 
Eq. ([2| that separates into rotation of the molecule as 
described by the 27r-periodic function £($) and torsion 
accounted for by the 27r-periodic function x(0d)- [Vig- 
orously, the bounds for depend on <pd [20 . At the 
time scales of interest, however, a molecule generally only 
makes a small fraction of a full rotation so that the effect 
of ^-<j)d coupling is negligible.] The separation is physi- 
cally motivated by considering the energy scales related 
to rotation and torsion. The energy scale of the prior is 
given by h 2 /(2I) = 1.3 jaeV. For torsion, on the other 
hand, the relevant energy is determined by the torsional 
potential, and a harmonic approximation of the poten- 
tial near the minimum at 39° yields a frequency corre- 
sponding to the energy 3.1 meV. Using that the period 
of motion is of the order Planck's constant divided by 
the energy, we therefore see that the molecule rotates (<£ 
change) on a nanosecond time scale, whereas the torsion 
(4>d change) is of picosecond duration. 

The small energy separation of the rotational levels 
compared to that of the torsional levels implies that many 
rotational states will be occupied in thermal equillibrium. 
Consequently, we expect the dynamics to behave classi- 
cally (see also Fig.[7|, and we only treat the torsional dy- 
namics quantum mechanically. The torsion is described 
by the stationary Schrodinger equation 
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where 0$, i = Br, F, is the angle of the i phenyl ring with 
respect to the kick pulse polarization axis, 1{ the mo- 
ment of inertia for rotation of the i phenyl ring around 
the stereogenic axis (J Br = 8911925, I F = 1864705) and 
Vt or ((j)Br — 0f) is the laser- free torsional potential as ob- 
tained from quantum chemical calculations (Fig. 3). By 
changing the coordinates to the dihedral angle (j)d = 
(j)Br — <j>¥ between the two phenyl rings and the weighted 
azimuthal angle $ = (1 — ?7)</>Br + ^f? characterizing 
the rotation of the molecule, with 77 = If /(If + ^Br), we 



with v = 1,2,... denoting the energy eigenstates. To 
solve this equation we expand the Hamiltonian onto an 
orthonormal basis of 27r-periodic functions, which we di- 
agonalize. The eigenstates below the torsional barriers 
exhibit an almost exact fourfold degeneracy, and linear 
combinations of the approximately degenerate states will 
also be stable at the time scales of the experiment. In 
particular, we can produce such stable states that are 
localized in the wells of the torsional potential. We will 
denote these states by |L^ in ), with i — 1, 2, 3, 4 and ^ m i n 
the smallest v of the four degenerate states. In Fig. [3] 
we show the first four energy eigenstates and the corre- 
sponding localized states along with the torsional poten- 
tial. The experiment is generally carried out on a gas of 
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FIG. 3: The field-free torsional states, (a) The four first (al- 
most) degenerate energy eigenstates lying 1.71 meV above 
the minimum of the torsional potential. From linear combi- 
nations of these we obtain the corresponding four localized 
states shown in (b) . The torsional potential is indicated with 
a dashed line. 



molecules in a mixed state given by the density matrix 



2=1 
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with the P's being weight factors that sum to unity, 
and where the latter sum includes the states above 
the torsional barrier. The kick pulse introduces dy- 
namics that we determine from the alignment poten- 
tial created by the pulse. By definition F(t) = Fo(t)X 
for the kick pulse, and we use a Gaussian envelope 
F (t) = F exp{-ln(2)[t/(r FW HM/2)] 2 / 2 } corresponding 
to an intensity full width at half maximum (FWHM) of 
tfwhm- We calculate the polarizability tensor in the MF 
frame as a function of the dihedral angle (see Table [i]) and 
consequently transform into the LF frame by the appli- 
cation of directional cosine matrices [21] . We then arrive 
at the interaction 

Vkick($,0d,*) = ~\f% (t)[a xx (^)cos 2 ($ + ^ d ) 

+ a yy (0 d ) sin 2 ($ + r}(j>d) 

- 2a xy (0 d ) cos($ + rj(f)d) sin(^> + rj(f) d )]. 

(5) 

For the succeeding analysis it is helpful to note a few 
qualitative features of the potential and we refer to Fig. [2] 
for graphical guidance. For a fixed dihedral angle the 
potential is minimal when the second most polarizable 
axis (SMPA), which is located 11° from the Br-phenyl 
ring and 28° from the F-phenyl ring, is parallel to the X 
axis. We will denote this by the || -geometry. Conversely, 
the potential is maximal if the molecule is rotated 90° 
from the || -geometry, and we will denote this by the _L- 
geometry. Next, in the || -geometry the potential favors 



a reduction of the dihedral angle, whereas an increase 
of the dihedral angle is resulting from the ^-geometry. 
Hence, the overall effect of the kick pulse will be to align 
the molecules into the || -geometry and drive a change of 
the dihedral angle. 

We proceed with the quantitative analysis of the field- 
induced dynamics. The state \xJ) is no longer an eigen- 
state of the torsion when the kick pulse is applied, but 
will develop into a wave packet 
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with the coefficients determined by 
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Once we find these new states of torsion the expecta- 
tion value of any operator O can be evaluated by tracing 
the product of the density matrix with the operator, i.e., 
(0(t)) — Tr[p(t)0]. In particular, we need the expecta- 
tion value of the dihedral angle and of the kick poten- 
tial from Eq. From the latter we obtain the torque, 
which causes the molecule to rotate into the || -geometry. 
If the molecule lies at an angle $ it is exposed to a torque 
— cKVkickt^, i))/d& directed along the Z axis and hence 
achieves an angular acceleration given by 



<9$ 



(8) 



Along with the initial conditions given at Eq. Q, the 
Eqs. ([7]) and (|sj) provide a set of coupled differential equa- 
tions that may be integrated to obtain the coordinates 
&(t) and (<pd(t)} at time t. Rather than solving these 
coupled equations, we, however, assume that the angle $ 
has the constant value ^>o during the short time interval 
of the kick pulse, and we integrate Eq. ([8| twice to arrive 
at 



$(£) = $ - 1 



-(- f 



dt'{V kick ($,t')) 



■ (9) 



Consistently, we solve Eq. ([7| with fixed at <3>o. 

For the present work we start out with a distribution 
of values for <l>o in the interval [— 7r, it] along with an ini- 
tial density matrix, p(to) and we propagate each member 
of the distribution according to Eqs. ^ and (|9|). This 
procedure yields an ensemble of $(t)'s and (^(t))'s and 
we finally employ the relations (/>Br(t) = $(£) + rj(<fid(t)} 
and 0f(£) = 3>(£) — (1 — rj)(<fid(t)} to find the angular 
distributions of bromines and flourines. 



B. Quantum chemical calculations 

To determine the field-free torsional potential we ob- 
tain the geometries of DFDBrBPh in various conforma- 
tions at the density functional theory level [22j [23] , with 
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TABLE I: The table lists the relevant polarizability components, ay , of DFDBrBPh in the MF frame as a function of the dihedral 
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B3LYP, Becke's three-parameter hybrid exchange func- 
tional [24] in connection with the Lee- Yang-Parr corre- 
lation [25] functional. We model the correlation of the 
uniform electron gas with the Vosko-Wilk-Nusair VWN5 
formulation [26], and we use doubly polarized triple-zeta 
quality basis-sets, TZVPP [27]. This level of theory has 
been shown to be suitable for the study of biphenyls [28] . 
To scan the potential energy surface of the relative tor- 
sional angle of the phenyl rings we optimize all coordi- 
nates except the torsional angle. 

We do all electronic structure calculations with the 
Turbomole 5.10 program suite [29]. In Table [j] we give 
the dynamic polarizabilities computed at a laser wave- 
length of 1064 nm. We note that the polarisabilities 
are quite stable with respect to choice of laser frequency. 
Using a value of 800 nm changes the individual compo- 
nents only minutely (maximally by 2%). The differences 
between static and dynamic polarisabilities are equally 
small. We calculate the vibrational frequencies analyt- 
ically [30], and a normal mode analysis, using internal 
coordinates, shows that the lowest frequency vibration 
corresponds almost purely to the torsion of the phenyl 
rings. We compute vibrational Raman cross-sections 
from derivatives of the polarizability tensor with respect 
to the normal modes of vibration. Table [ill shows the 
frequencies and Raman cross sections of the lowest vi- 
brations of DFDBrBPh, computed at B3LYP/TZVPP 
level. The results provide the justification for assuming 
that the kick pulse primarily interacts with the lowest 
normal mode, since its Raman scattering cross section 
is over four times higher than the second lowest normal 
mode and 20 times higher than any other normal mode. 
We note that the character of the second lowest (and sec- 
ond strongest Raman active) mode contains no torsional 
motion. 



TABLE II: Vibrational frequencies and relative Raman cross 
sections of the eight lowest lying normal modes of DFD- 
BrBPh, computed at B3LYP/TZVPP level. 



mode 


frequency (cm 1 ) 


cross section (%) 


v\ 


23.9 


100.0 


V2 


55.7 


21.4 


V3 


60.7 


2.8 


l/ 4 


111.9 


5.1 


"5 


145.2 


1.3 


vq 


194.7 


3.6 


h> 7 


195.3 


1.3 


V8 


213.4 
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Na2C03 (3 mL, degassed). We next heat the mixture 
to 90° C under argon for 5.5 h and cool it to room tem- 
perature. The mixture is diluted with water and ex- 
tracted three times with Et2 0. We dry the combined 
organic fractions with MgSC^, filter and concentrate to 
give approximately 1.4 g of a crude mixture. Based on 
NMR we find that this mixture contains both the de- 
sired mono-Suzuki coupling product [31], double-coupled 
products, and a homo-coupling product from the boronic 
acid. Careful column chromatography on Si02 eluting 
with n-hexane and c-hexane finally allows us to isolate 
182 mg 3,5-dibromo-3',5'-difIuorobiphenyl (DFDBrBPh) 
as a white amorphous solid (contains <3% of the homo- 
coupling product DFDBrBPh). NMR spectroscopic data 
(see also Fig.|§: X H NMR 5 (400 MHz, CDC1 3 ) 7.69 (t, J 
1.7 Hz, 1H), 7.61 (d, J 1.7 Hz, 2H), 7.05 (m, 2H), 6.85 (m, 
1H). 13 C NMR S (100 MHz, CDC1 3 ) 163.5 (2C), 142.3, 
141.5, 133.8, 128.9 (2C), 123.5 (2C), 110.1 (2C), 103.8. 
Elemental analysis, calculated: C 41.42%, H 1.74%, Br 
45.92%; found: C 41.20%, 1.66%, Br 45.90%. 



IV. EXPERIMENTAL METHODS 

A. Synthesis of 3,5-dibromo-3',5'-difluorobiphenyl 

To a round-bottom flask equipped with a magnetic 
stirring bar under argon we add 3,5-dinuorophenyl- 
boronic acid (474 mg, 3.0 mmol), 1,3,5-tribromobenzene 
(1134 mg, 3.6 mmol), and Pd(PPh 3 ) 4 (69.6 mg, 2 mol%). 
The flask is then evacuated and backfilled with argon 
twice. To this mixture of solids we add toluene (18 
mL, degassed), EtOH (6 mL, degassed), and 2M aq. 



B. Setup of the laser experiment 

A schematic of the experimental setup is shown in 
Fig. [5] Less than hundred milligram of solid DFDBrBPh 
is heated to 130° C and expanded in 90 bar of He into 
vacuum using a pulsed Even-Lavie valve [19] [32]. The 
resulting supersonic molecular beam is skimmed and sub- 
sequently crossed by three, focused, pulsed laser beams at 
90°. The first laser beam originates from a Q-switched 
Nd:YAG laser and consists of 9-ns-long pulses at 1064 
nm. It is focused to a spotsize coq ~ 35 jam yielding a 
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FIG. 4: (Color online) (a), Synthesis of DFDBrBPh by Pd- 
catalyzed cross-coupling, (b)-(d), NMR-spectra of synthetic 
DFDBrBPh. Integrated signals are indicated by green (grey) 
curves, and the corresponding numbers are situated directly 
below the respective curves. 



peak intensity of ~ 7x 10 11 W/cm 2 . The purpose of these 
pulses is to align the stereogenic axis of the DFDBrBPh 
molecules. The second laser beam consists of 0.7-ps-long 
kick pulses (A = 800 nm), obtained by passing part of 
the output from an amplified Ti: Sapphire fs laser system 
through a grating stretcher. The spotsize of c^o = 43 
/im yields a peak intensity of ~ 5 x 10 12 W/cm 2 . These 
kick pulses are used to induce torsion and rotation of the 
molecules. The third laser beam consists of 25-fs-long 
pulses, which we achieve by passing another part of the 
output from the Ti: Sapphire system through a hollow 
wave guide compressor setup. The pulses have a spotsize 
ujo = 25 /im and a peak intensity of ~ 2 x 10 14 W/cm 2 . 
These probe pulses Coulomb explode the molecules and 
the produced ions are extracted by a weak static electric 
field, in a velocity imaging geometry, and projected onto 
a two-dimensional detector, consisting of a micro channel 
plate (MCP) detector backed by a phosphor screen. We 
detect the F + and Br + ions separately by time gating 
of the MCP detector around the respective arrival times 
of the ions. We record the ion images on the phosphor 
screen by a CCD camera and determine the coordinates 
of each individual ion hit. Figure [6ja) displays the ion 
images of F + and Br + at selected probe times. 
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FIG. 5: (Color online) Schematic representation of the ex- 
perimental setup, showing the vacuum system and the three 
laser beam paths. The Nd:YAG (alignment) pulses and the 
probe pulses are polarized vertically, i.e., perpendicular to the 
detector plane, whereas the fs kick pulse is polarized in the 
plane of the detector. 
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FIG. 6: (Color online) (a) Ion images of F + and Br + frag- 
ments at probe times t p . The ns pulse is polarized per- 
pendicularly to the image (detector) plane and the 5 x 10 12 
W/cm 2 , 0.7 ps (FWHM) kick pulse is polarized horizontally, 
(b) Angular distribution of the F + ions, at t p =1.47 ps, ob- 
tained by radially integrating the corresponding F + ion im- 
age. The splitting of the pairwise peaks is twice the av- 
erage angle, (<Pf+- kick. P oi)> between the F + ion recoil (and 
thus the F-phenyl plane) and the kick pulse polarization, (c) 
(0f+- kick. P oi) as a function of t p , for times where a clear four- 
peak structure is visible in the angular distributions. The 
curve is a fit of the sum of a linear and a harmonic function 
to the experimental points (squares). 



V. RESULTS AND DISCUSSION 

First, we establish that the ns pulse is able to hold 
the C-C bond axis of DFDBrBPh along its polarization 
direction. At t p = —0.87 ps the F + image, displayed in 
Fig. [6^a), is almost circularly symmetric and the small 
deviation from circular symmetry is explained by noting 
that the kick pulse has a finite value at —0.87 ps. The 



absence of ions in the innermost region is only compati- 
ble with the C-C bond axis being aligned perpendicular 
to the detector plane and the F-phenyl planes uniformly 
distributed around the C-C bond axis. The correspond- 
ing Br + image is also essentially circularly symmetric al- 
though the hollow-out of the center is not as pronounced 
as for F + . 



7 



Next, we investigate the effect of the kick pulse. Al- 
ready at t p = 0.47 ps the circular symmetry is broken 
and both F + and Br + ions start to localize around the 
polarization direction of the kick pulse. The F + ions re- 
main within the radial range of the ring structure, which 
shows that the kick pulse does not cause any significant 
distortion of the C-C bond axis alignment but rather ini- 
tiates an overall rotation of the molecule around this axis. 
The latter is expected since the torque imparted by the 
kick pulse will force the SMPA (Fig. 2) to align along the 
kick pulse polarization after a delay determined by the 
kick strength. At t p = 1.47 ps the localization sharpens 
for both ions species, but whereas the Br + ions are sim- 
ply more localized around the polarization axis the F + 
ion distribution exhibits a four dot structure. This be- 
haviour is compatible with alignment of the SMPA along 
the kick pulse polarization. In the limit of perfect SMPA 
alignment, the Br + ions would appear as two pairs of 
strongly localized regions on the detector corresponding 
to molecules with their Br-phenyl plane located either 
11° clockwise or counterclockwise to the SMPA. Each 
pair of Br + ions would then be angularly separated by 
22°. In practise and consistent with theory (see Fig. [8| 
the SMPA alignment is not strong enough to resolve the 
two Br + regions. It is, however, sharp enough to resolve 
the two pairs of F + ion hits due to the much larger offset 
(28°) of the F-phenyl plane from the SMPA. At t p = 2.47 
ps the Br + distribution has localized further showing that 
the Br-phenyl plane has rotated into stronger alignment 
with the kick pulse polarization. If the torsion was un- 
affected, i.e., (j)d remained unchanged, the F + ion image 
at t p = 2.47 ps should exhibit a distinct four-dot struc- 
ture similar to the image at t p = 1.47 ps but with a 
larger angular splitting between the ion regions in each 
of the two pairs. Clearly, the four-dot structure at 2.47 ps 
is significantly blurred compared to the case at 1.47 ps. 
Thus, we conclude that the kick pulse not only sets the 
molecule into controlled rotation around the C-C axis, it 
also initiates torsional motion. 

Beyond 3.8 ps the angular confinement of both the 
F + and the Br + ions is gradually lost due to continued 
rotation around the C-C axis with dihedral dynamics im- 
posed, and the images for both ion species eventually re- 
gain their circularly symmetric form (see images at 12.15 
ps in Fig. [6^a)) identical to the pre-kick pulse images at 
-0.87 ps. 

Further insight into the effect of the kick pulse is ob- 
tained by analyzing the angular distribution of the F + 
ions as a function of t p . An example of the F + angular 
distribution at t p = 1.47 ps is shown in Fig. [6^b). Fig- 
ure [6jc) displays the average angle between the F-phenyl 
rings and the kick pulse polarization as a function of t p 
in the time interval where a clear four-peak structure is 
visible in the angular distributions. The increase from 
22.5° at 1.47 ps to 26.5° at 3.8 ps shows that the F- 
phenyl plane gradually moves away from the kick pulse 
polarization due to the overall rotation of the molecule 
around the C-C bond axis, and we ascribe the recurrent 
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FIG. 7: (Color online) Comparison of the classical model for 
rotation at various intensities with the experiment of Ref . [33] . 
The theory captures the laser induced rotation within the 
first 12 ps, but generally overestimates the degree of angular 
confinement. 



dips to a periodical variation in <\>&. We estimate the pe- 
riod to be ~ 1 ps and the amplitude to ~ 0.6° for this 
oscillation. 

Before presenting the theoretical results on DFD- 
BrBPh, we want to confirm the validity of the classical 
treatment of rotation. To this end we apply the model 
to a test case, where the simpler 3,5-diflouroiodobenzene 
(DFIB) has its most polarizable axis fixed by a ns laser 
pulse and is then set into rotation about this axis by an 
orthogonally polarized femtosecond pulse. Here, there 
is no torsion and the rotational dynamics have been ob- 
served directly recently [33]. In Fig. [7| we compare the 
theoretical and experimental results and find good agree- 
ment. 

Figure [8] shows the results of a calculation with laser 
parameters identical to the experimental values and an 
initial rotational temperature of K. Prior to the kick 
pulse the angular distributions of the Br- and F-phenyl 
rings are isotropic as in the experiment. Maximum align- 
ment of the SMPA occurs at 1.3 ps and the confine- 
ment of the F-phenyl rings at a large angle with re- 
spect to the kick pulse polarization (cf. middle panel, 
Fig. |8ja)) explains the distinct four-dot structure ob- 
served at tp = 1.47 ps in the experimental F + ion image 
(Fig. [6^ a)). Also, at t p = 1.47 ps the confinement of the 
Br-phenyl rings at a small angle with respect to the kick 
pulse polarization predicts a much less distinct, if any, 
four-dot structure in good agreement with the Br + ion 
image. At t p = 2.47 ps the angular localization of the F- 
phenyl rings has broadened (right panel, Fig.[8ja)) and a 
blurred four-dot structure is seen, consistent with the ex- 
perimental result at t p = 2.47 ps. The distribution of the 
Br-phenyl rings is also broadened (right panel, Fig.^b)), 
but remains localized around the kick pulse polarization 
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FIG. 8: (Color online) Angular distributions of (a) F-phenyl 
and (b) Br-phenyl rings at t p = 0.08, 1.47 and 2.47 ps. (c) 
Expectation value of the dihedral angle for a molecule starting 
out with the SMPA aligned along the kick pulse polarization. 
The kick pulse is as in Fig. [6] 



fully consistent with the Br + ion distribution, recorded 
at 2.47 ps. 

The theoretical value (0^) exhibits oscillations 
(Fig. |8][c)) with a period of ~ 1.2 ps and amplitude of 
~ 2.45°. The period agrees well with the experimental 
value (~ 1 ps), and smaller modulation in (<^) is ex- 
pected in the experiment (~ 0.6°) partly since here the 
SMPA is not pre-aligned and partly due to deviations 
from the stated laser peak intensity. The behavior is as- 
cribed to a wave packet of normal modes in the torsional 
double well potential (Fig. [3^b)) for a molecule starting 
out with the SMPA aligned. The qualitative agreement 
of Figs. |6|b) and [8jc) corroborates the interpretation of 
the kick pulse inducing time-dependent torsional motion. 
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FIG. 9: (Color online) Time evolution of the dihedral angle 
for a molecule starting out as (a) an R a or (b) an S a enan- 
tiomer. Initially, the molecule is 3D oriented with the Br- 
phenyl end pointing out of the paper and the SMPA aligned 
at an angle of 13° with respect to the kick pulse polariza- 
tion. The kick pulse triggering the torsional motion has a 
peak intensity of 1.2 x 10 13 W/cm 2 and duration (FWHM) of 
1.0 ps. The torsional potential is scaled down by 1/4 rather 
than by increasing the kick strength. In practise reduction 
of the torsional barrier is possible, for instance, by modifying 
the aromatic rings or by using halogen substituted bipheny- 
lacetylene. 



VI. PERSPECTIVES 



We will now discuss some perspectives of the demon- 
strated laser control of torsion. For one thing, we sug- 
gest a time-resolved study of de-racemization [34, 35j[36j 
[37| 138] , where one enantiomer is steered into its mirror 
form. To do this we break the inversion symmetry of the 
C-C bond axis by orienting each molecule. This allows 
us to discriminate between the two enantiomeric forms. 
Theory [39 , 40] and experiment [JTJ|42] show that orien- 
tation can be added to 3D alignment by combining the 
ns alignment pulse with a static electric field. Next, the 
excitation of torsional motion is optimized by eliminating 
initial rotation around the C-C bond through alignment 
of the SMPA (Fig. 2) prior to the kick pulse using an el- 
liptically rather than a linearly polarized ns pulse [7[ |43] . 
Finally, the interaction strength between the molecule 



and the kick pulse needs to be increased either through 
higher intensity, a longer kick pulse or by trains of syn- 
chronized kick pulses [44] [45j [46] . 

Assuming initial orientation and confinement of the 
SMPA, we have calculated the <pd dynamics for both con- 
formations the DFDBrBPh molecule with a reduced tor- 
sional barrier (see caption to Fig. [9]). The results are 
shown in Fig. [9] and clearly, the present method would 
allow for a time-resolved study of the transition from one 
enantiomer into the other. A quantitative analysis of the 
efficiency of the process shows that after the pulse 99% 
of the molecules starting out as R a has changed into S a 
enantiomers, whereas only 13% of the S a enantiomers 
changed into R a . The inverse process causing an excess 
of R a enantiomers, is achieved simply by inverting the 
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orientation of the molecules. 

We believe that the results presented in this paper also 
fuel the field of molecular junctions [47, 48, 49 with new 
possibilities. The conductivity of molecules like biphenyl 
will depend on the dihedral angle [50 J, and as such this 
type of molecules, if placed between two leads, can func- 
tion as a molecular junction [5TJ [52j [53j [54] that may 
be used to control the charge flow and act as, e.g., a 
switch [43] . Importantly, we can realize quantum control 
of such a system by means of kick pulse schemes and in 
this way tailor very specific torsional wave packets [55] 
to dictate the detailed time-dependent dihedral motion 
and thereby the current flow through the molecular junc- 
tion. Seen through the light of controllability, molecular 
junctions based on laser controlled torsion complements 
the available schemes, such as the mechanical break junc- 
tions [56 and the previously suggested resonantly light 
driven molecular junctions [57] [58j [59] . 

VII. CONCLUSION 

We presented a fs time-resolved study on laser con- 
trolled torsion of axially chiral molecules. Experimen- 
tally, the symmetry axis of DFDBrBPh molecules was 
held fixed in space by a long alignment laser pulse while 



a fs kick pulse, polarized perpendicular to the fixed axis, 
was applied. The torsion as well as the overall rotation of 
the molecule was monitored with fs time resolution using 
a delayed probe pulse. To substantiate the observations, 
we developed an original theoretical model, which pro- 
duces results in good agreement with the experiment. 
The model divides the dynamics of the molecule into a 
the rotation with classical behavior and quantum me- 
chanical torsional dynamics. As such the model provides 
a transparent physical intepretation of essential features 
of laser controlled torsion of DFDBrBPh. Based on the 
theoretical model, we further discussed the extension of 
the experiment to a setup that will facilitate a fs time- 
resolved study of de-racemization. Finally, we pointed to 
perspectives of using molecules like DFDBrBPh in molec- 
ular junctions in order to realize laser controlled charge 
flow. 
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